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INTRODUCTION 


To complete the requirements of contract NAS8-31904, which 
was initiated to investigate lubrication effects on bearing thermal 
performance, an investigation was performed to determine the feasi- 
bility of using the SKF program SHABERTH for simulating the perform- 
ance of cryogenically lubricated ball bearings. As a part of this 
study, the particular application chosen for SHABERTH was to simu- 
late the performance of the Space Shuttle main engine turbo-pump 
and pre-burner bearing system. 

A concurrent effort was made by cognizant MSFC computation lab 
personnel to convert SHABERTH to accept standard engineering English 
units, but this effort has not been completed at this writing. The 
units used in SHABERTH are the international system of units. 
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TECHNICAL DISCUSSION 


The particular problems that SHABERTH has been designed to 
solve are: 

• Elastic contact deformations 

• Diametral clearance change 

• Heat generation temperature distribution and 
friction effects with lubrication 

According to the programmers/users manual (Reference 1) , the 
program is applicable to bearing systems with up to 5 roller and/ 
or ball bearings. The application chosen to evaluate the SHABERTH 
program was the Space Shuttle main engine turbopump and pre-burner 
bearings. A schematic of the bearing system is shown in Figure 1. 
All four bearings are angular contact ball bearings, which are 
lubricated by liquid oxygen and are subjected to both axial and 
radial loads. The bearings are spring loaded one against the other. 

A simplified analytical model was initially developed and was 
subjected to several loading cases. The loading conditions and 
results output by SHABERTH are shown in Table 1. Although the 
program did not indicate an error might exist, some of the results 
appear to be unrealistic. A more detailed model was developed, 
which included the lubrication effects of liquid oxygen. Attempts 
to execute this model always ended in an error return from a SHA- 
BERTH subroutine. Telephone contact with SKF programmers, confirmed 
the observation of existing errors in the thermal portion of the 
program as well as program code which make it impossible to use 
cryogenic temperatures in the model. During the course of this 
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study, after scores of computer runs, it was found that SHABERTH 
could not successfully solve this problem when both radial and 
axial loads are applied. It was also ascertained by trial and 
confirmed by SKF programmers that SHABERTH could not simulate a 
bearing system in which the bearings are pre-loaded one against 
the other, because axial loads must be applied to the shaft, and 
they are algebraically summed. In this case, the shaft thrust 


load would be zero. 



MODEL DEVELOPMENT 


Although the user's manual gives reasonably clear instructions 
on input data development, there are a few areas that need to be 
expanded upon. In order to do a bearing diametral clearance change 
analysis, the thermal expansion is referenced to cold conditions 
at 68 °F. The input parameter for contact angle is not the manu- 
facturer's specified contact angle but is an auxiliary contact 
angle. This angle is determined by the following equations: 


(1) 

( 2 ) 

(3) 

(4) 

(5) 


t 

A 


= COS 


-1 


2A-Pd 

2A 


= r o + r i 


D 


Sd = [Pd - 2A(1 - coscA 0 ) ]/cos^ 0 
by rearranging terms 


Sd / cos 2 (X Q - A cos<^ 0 + A - Pd/2 = 0 

then by substitution 

Pd - 2A (1-cos c( ) Pd 

2_ - a cos''/ q + A - — 


0 


where , 

oC ~ manufacturer's designated contact angle 
r Q = outer raceway groove radius 
r^ = inner raceway groove radius 
D * ball diameter 
Pd = diametral clearance 
do - input parameter for SHABERTH 
(solved by iteration) 

Sd = auxiliary clearance parameter 
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The auxiliary contact angle (U Q ) should be less than the manu- 
facturer ' s contact angle . With these introductory remarks about 
model development, the following statements describe additional 
idiosyncrasies and/or important or unclear input parameters en- 
countered in a step by step model development. 

A SHABERTH model consists of data organized in four distinct 
categories. The categories are title cards, bearing data cards, 
thermal data cards, and shaft data cards, respectively. The "Title 
Data Group," consisting of only two cards, is clearly defined, 
except that data should either have a decimal value or be right 
justified in the field, as is true for all SHABERTH data. 

The "Bearing Data Group" consists of sets of up to 16 cards 
for each bearing in the model. The sets must be input sequentially 
in the order of increasing distance from the origin along the "X" 
axis. Card type 2 contains the input parameter discussed earlier, 
the auxiliary contact angle. Card types 3, 5 and 6 are not input 
for ball bearings. Card types 7 and 8 are not input for NPASS = 0, 
and card types 9 through 14 are deleted from the input set if the 
diametral clearance change analysis is not to be performed. If the 
clearance change is to be performed, the parameters for inner and 
outer ring mean diameter must be input. This is an interesting 
input parameter since the inner ring has only one shoulder. Figure 
2 demonstrates the determination of these parameters. Card types 
15 and 16 are deleted from the input set for NPASS = 0 . It should 
be noted that card types 15 and 16 are for lubricant properties, 
and it has already been stated that SHABERTH would not run with the 
cryogenic thermo-physical properties of liquid oxygen. 



The "Thermal Data Group" consists of nine card types, with 
some types repeated for each bearing. Figure 3 provides the node 
numbering used for the thermal model developed. SHABERTH also 
errored in every attempt to execute the thermal model. Cognizant 
SKF programmers were contacted and some updates to correct some 
of the problems were received and have been given to cognizant 
MSFC computation lab support personnel to incorporate into the 
base line program to be maintained by the computation lab. The 
input for card types 7 and 8 are confusing du<_ to the terminology 
used in the user's manual. What is intended is that you can define 
10 solid conductance values (indicated by identification numbers 
1 through 10) , 10 free convection conductance factors (indicated by 
identification numbers 11 through 20), et cetera, see Appendix I. 

The •'Shaft Data Group" consists of three card types which 
describe the shaft geometry, bearing locations on the shaft and 
shaft loading respectively. Card type 1 is used to describe the 
shaft wall thickness changes, with up to 20 changes in diameter 
allowed (20 cards) . Card type 2 locates the bearings along the 
X axis of the shaft, 1 card per bearing. There are two types of 
"type 3" input cards, all shaft loads in the X-Y plane are input 
and terminated by a blank card, then all shaft loads in the X-Z 
plane, followed by a blank card. Even if there is no loading in 
the X-Y plane, or X-Z plane, the type 3 card must be input, fol- 
lowed by a blank card. 
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RESULTS AND CONCLUSIONS 


The computer generated results for the thermal model are 
contained in Appendix I. Several attempts were made to get the 
model to run, but as previously stated, errors in the program 
always caused the execution to terminate. Hopefully the cor- 
rections sent by SKF, which have been given to cognizant MSFC 
computation lab contractor personnel, will eliminate the error 
termination. However, the user should remember that the program 
bases the diametral clearance change analysis on an initial 
temperature of 68°F. Also the program expects the lubricant 
viscosity data to be input at 100°F and 210°F. 

Appendix II contains the computer generated output for the 
diametral clearance change model. According to the user's manual, 
a positive contact angle allows the bearing to accept a positively 
directed axial load transmitted by the shaft. However, I changed 
the contact angles to negative on the second and fourth bearings 
on a recommendation from SKF programmers, because the program 
could not find a solution when both radial and axial loads were 
applied to the model. I overrode the default maximum iterations 
and allowed up to 1000 iterations, with no solution. SKF program- 
mers stated that they thought that they had programmed SHABERTH 
to avoid getting caught in infinite iterations, but they obviously 
have not. The runs always terminated because of maximum CPU time 
of 30 minutes. When the radial load was eliminated, the program 
converged in only 7 CPU seconds. The applied loads were used only 
as part of a parametric study, and not as a realistic loading case. 
It is obvious, however, that the fourth bearing does carry part 
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of the load even though it has a negative contact angle. The 
user's manual states that the "solution of a single# radially 
and axially loaded problem is impossible." Another impossible 
situation described in the user's manual is a problem in which 
a radial load and an outer ring misalignment and a zero applied 
moment are specified. It does not imply that there may be other 
combinations of problems which the program cannot solve. Obviously 
there are such combinations as described above. 

Although numerous problems have been encountered in the utili- 
zation of SHABERTH, it is obvious from the complexity of the pro- 
gram and the limited experience gained during this study that 
SHABERTH has the potential for solving certain bearing problems. 

It is recommended that a continuing effort be expended to perform 
a series of parametric studies of various loading conditions, 
utilizing the existing shuttle engine model. 
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Figure 1. Space Shuttle Pre-burner and Turbopump Model Schematic 
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Outer Mean Diameter * 13.420 + 3.19 + (2) (.265) ]/2. * 3.57 
Inner Mean Diameter * (2.694 + 2.904J/2. 

(See Reference 2.) 

Figure 2. Determination of Mean Diameters 
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NOTES: 


Each bearing was incremented by (5), i.e., the second 
bearing node numbers were 6 through 10, the third, 11 
through 15, and the fourth, ]6 through 20. The liquid 
oxygen was node 25. 1 is the inner ring, 2 are the 

balls, 3 is the cage, 4 is the outer ring, 5 is *he 
shaft, and 21 through 24 are housing nodes. 


Figure 3. Thermal Model Node Designations 
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CASE #1 


Parameter 


Bearing Bearing Bearing Bearing 


Shaft 


Contact Angle 


Radial Load 


Thrust Load 


Lubricated 


Thermal 


Clearance 

Change 


Total Thrust Load 


T initial 0< - 


T f inal 


22.5 


1779.3 



Cage Speed RPM 


Max Orbit RPS 


Max Outer 
Contact Angle 


Max Inner 
Contact Angle 


Max H z Stress 
N/MM 2 

Outer 


Inner 


2093 


Table 1. Initial SHABERTH Results 


O Indicates Ball Number 







































































Bearing 

Bearing 

Bearing 

Bearing 

Parameter 

# 

it 

# 

« 


1 

2 

3 

4 

Contact Angle 

24.34 

24.34 

20.5 

20.5 

Radial Load 

1779.3 

1779.3 

5337.8 

i 

5337.8 


Shaft 


Thrust Load 


Lubricated 


Thermal 


Clearance 

Change 


Total Thrust Load 


C 


T final ° c 


S R -F. 


D x 


D. 


D, 



Life Hours 
Cage Speed RPM 


Max Orbit RPS 


Max Outer 
Contact Angle 


Max Inner 
Contact Angle 


Max Hz Stress 
N/MM 2 


Outer 


Inner 



1766 


-.253 


.129 


1 x 10" 9 


1 x 10 -5 


-257. 


491 


11610 


1216 ® 


.19 


.29 


2134 w 


2300 ® 


-264. 

460 

11610 


1216 ® 


.19 


.29 


320. 

322. 

68.2 

64.9 

11810 

11810 

1236 ® 

1236 ® 

-.07 

1 

• 

o 

-^1 

-.08 

-.08 



2142 


2764 


2764 


2562 w 


2779 ^ 


Table 1. (Continued) 



































































Parameter 


Bearing 

# 

1 


Contact Angle 

24.34 

Radial Load 

1779.3 

Thrust Load 

1890.45 

Lubricated 

No 

Thermal 

No 

Clearance 

Chanqe 

No 

Total Thrust Load 


T initial ° c 

-176. 

T final °C 

-176. 

Z R.F. 

2092 

D x 

-1 x 10" 3 

D y 

7 x 10~ 2 

d z 

5 x 10 -9 

M y 

-1 x 10" 3 

M Z 

-22250 

Life Hours 

1803 

Cage Speed RPM 

12800 

Max Orbit RPS 

mm 

Max Outer 
.Contact Anqle 

17.96 ® 

Max Inner 
Contact Anqle 

45.03 ® 

Max Hz Stress 
N/MM 2 


Outer 

1828^) 

Inner 

1894 OD 






































































CASE #4 


Parameter 

Bearing 

# 

1 

.. _ . 

Bearing 

# 

2 

Bearing 

# 

3 

Bearing 

« 

4 

Shaft 

Contact Angle 

24.34 

24.34 

20.5 

20.5 


Radial Load 

1779.3 

1779.3 

5337.8 

5337.8 


Thrust Load 

3780.9 

3780.9 

0 

0 


Lubricated 

No 





iBBBBMBB 

hhhhhi 

B BB 

Thermal 

No 





B 

Bi B 

B IBM 

Clearance 

Chanqe 

No 

r 




B B 

B 

B Bfli 

Total Thrust Load 





7562 


Bn 

bbi b 

BB 

^initial °C 

-176. 





T final ° c 

-176. 





BggjjpBBili 




En.r. 

2208 

2378 

7750 

9464 


D X 

2.3xl0” 2 








D y 

5 . 9xl0~ 2 

6 . OxlO -2 

8.1x10 2 

8.8x10-2 


D Z 

msM 





■1 

■IB 

B£in$f$£ii£BB 

My 

-lxlO" 3 

-1x10-3 

-5x10-3 

-6xl0~3 


M Z 

-24060 

-25770 

-89180 

-106600 


Life Hours 

2109 

1722 

88.4 

B 


Cage Speed RPM 

13810 

13810 

12510 

12470 


Max Orbit RPS 

RR^ 

|RnB[ 

Qyyi 

1408 @ 


Max ''uter 
.Cont.uct Angle 

19.44 © 

19.90 © 

17.56 © 

17.68 © 


Max Inner 
1 Contact Anqle 

| 

56.92 © 

39.12 ® 

38.94 © 


Max H z Stress 
N/MM 2 






Outer 

1779 © 

1808 © 


2553 © 


Inner 

1825 © 

1881 © 

2598 © ; 

2777 © 

i 


Table 1. (Continued) 
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CASE *5 



Bearing 

Bearing 

Bearing 

Bearing 


Parameter 

# 

# 

# 

# 

Shaft 


1 

2 

3 

4 


Contact Angle 

24.34 

24.34 

20.5 

20.5 


Rar ial Load 

1779.3 

1779.3 

5337.8 

5337.8 


Thrust Load 



2668.9 

2668.9 


Lubricated 

No 









Thermal 

NO 









Cleara) ce 





5338 

Chanqe 






Total Thrust Load 

-176 °C 










-176°C 









T final ° c 

-176 °C 










2 R - F * 

2017 

2158 

7091 

8308 


D v 

-3.X10" 2 





X 





_ D -y __ 

8 x 10 -2 

8 x 10" 2 

.101 

.105 


d 2 

4 x 10” 9 

4 x 10" 9 

6 x 10" 9 

6 x 10~ 9 


M y 

-9xl0“ 4 

-lxlO -3 

p4 x 10“ 3 

-4 x 10 -3 


M z 

-19540 

-20690 

-69150 

-79560 


Life Hours 

1411 

1185 

78.6 

50.5 


Cage Speed RPM 

12250 

; — 

12230 

| . . 

12110 

(- -- 

12070 


Max Orbit RPS 

1 

1341 ® 

1339 ® 

1296 ® 

1289 ® 



15.97 ® 

16.24 ® 

14.02 ® 

14. 07*® 


Max Inner 
Contact Anqle 

34.33 ® 

34.15 O) 

26.20 ® 

25.37 ® 


Max H z Stress 
N/MM 2 






Outer 

1893 ® 

1919 ® 

2472 ® 

2577 ® 


Inner 

1982 CO 

2029 ® 

2670 ® 

2805 ® 



table 1 . (Continued) 
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APPENDIX I 


THERMAL MODEL WITH LOX 


THERMO-PHYSICAL PROPERTIES 
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